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This paper exploits the fact that Hamiltonian flows associated with a time-independent 
H can be viewed as geodesic flows in a curved manifold, so that the problem of stability 
and the onset of chaos hinge on properties of the curvature K a f, entering into the Jacobi 
equation. Attention focuses on ensembles of orbit segments evolved in representative two- 
dimensional potentials, examining how such properties as orbit type, values of short time 
f-» Lyapunov exponents x, complexities of Fourier spectra, and locations of initial conditions 

on a surface of section correlate with the mean value and dispersion, (K) and ag, of the 
' (suitably rescaled) trace of K a b. Most analyses of chaos in this context have explored the 

effects of negative curvature, which implies a divergence of nearby trajectories. The aim 



here is to exploit instead a point stressed recently by Pettini (Phys. Rev. E 47, 828 [1993]), 
namely that geodesies can be chaotic even if K is everywhere positive, chaos in this case 
arising as a parameteric instability triggered by regular variations in K along the orbit. For 
ensembles of fixed energy, containing both regular and chaotic segments, simple patterns 
exist connecting (K) and for different segments both with each other and with the short 
time x : Often, but not always, there is a nearly one-to-one correlation between (K) and 
o~g, a plot of these two quantities approximating a simple curve. Overall x varies smoothly 
along the curve, certain regions corresponding to regular and "confined" chaotic orbits 
where x ls especially small. Chaotic segments located furthest from the regular regions 
tend systematically to have the largest x's. The values of (K) and (and in some cases 
x) for regular orbits also vary smoothly as a function of the "distance" from the chaotic 
phase space regions, as probed, e.g., by the location of the initial condition on a surface 
of section. Many of these observed properties can be understood qualitatively in terms of 
a one-dimensional Mathieu equation, in which parametric instability is introduced in the 
simplest possible way. 
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I. INTRODUCTION AND MOTIVATION 



It is well known [1] that the flow associated with a time-independent Hamiltonian H = 
\5 ab PaPb + V{x a ) can be reformulated as a geodesic flow in a curved, but conformally flat, 
manifold. Specifically, let 

ds 2 = W(x a )6 ab dx a dx a , (1) 

where E is the conserved energy associated with the time-independent H and the conformal 
factor 

W{x a ) = E- V{x a ) (2) 

is equal numerically to the kinetic energy associated with a trajectory at the point x a . It 
then follows that, with the further identification 

ds = V2Wdt, (3) 

the geodesic equation for motion in the metric g a b = W5 a b is completely equivalent to the 
Hamilton equations 

dx a _ dH ^ dpa _ _dH_ 

This implies that the confluence or diverence of nearby trajectories x a {s) and [x + ^] a (s) 
is determined by the Jacobi equation, i.e., the equation of geodesic deviation, which takes 
the form 

^ = -fi; c /A^-r c f, (5) 

where R a bcd is the Riemann tensor associated with g a b and D/Ds = u a V a denotes a direc- 
tional derivative along u a = dx a /ds. Linear stability or lack thereof for the trajectory x a (s) 
is thus related to R a bcd or ) more precisely, to the curvature K a c . If, e.g., R a bcd is everywhere 
negative, so that K a c always has one or more negative eigenvalues, the trajectory must be 
linearly unstable. 

It would seem intuitive that, if the curvature of g a b is everywhere negative, so that nearby 
trajectories always tend to diverge, every geodesic will behave in a fashion that is manifestly 
chaotic. If one assumes that the manifold is compact, so that trajectories are restricted to 
a region of finite volume, this intuition can be elevated to a theorem. For example, geodesic 
flows on a compact manifold with constant negative curvature are necessarily iT-flows, 
where generic ensembles of initial conditions evolve towards a microcanonical distribution 
at a rate set by the magnitude of the curvature [2] . If the curvature is everywhere negative 
but not constant, the flow is more complex, but one can still infer [3] a chaotic evolution 
towards a microcanonical distribution at a rate bounded from below by the least negative 
value of the curvature. 

When the curvature is not everywhere negative, much less is known. Nevertheless, the 
preceeding has motivated the expectation that, in many dynamical systems, chaos should 
be associated with (regions of) negative curvature. In particular, several authors (cf. [4,5]) 
have sought to use negative curvature to explain the fact that the gravitational iV-body 
problem for a large number of objects of comparable mass is chaotic in the sense that the 
evolution manifests an exponentially sensitive dependence on initial conditions. 

However, as stressed recently by Pettini [6], not all chaos can be associated with negative 
curvature. In particular, one can have large measures of chaotic orbits even for systems and 
energies where K a b is everywhere positive. In retrospect, this is easy to understand. Viewing 
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the Jacobi equation as a matrix equation, one can solve at any given point in space to derive 
eigenvectors {Xj} and eigenvalues {Aj}, each pair solving a linear equation of the form 

= -A.X,. (6) 

When the curvature is everywhere positive, \ > 0, so that the solutions are oscillatory 
rather than exponential. If the Aj's were constant along the trajectory, one could thus infer 
stable oscillations. In general, however, the Aj's are not constant, depending instead on the 
unperturbed x l {s) since R a bcd an d u a both change along the trajectory. It follows that, even 
assuming an everywhere positive curvature, eq. (6) should be interpreted as an oscillator 
equation 

^ = (7) 

The obvious point is that solutions to eq. (7) can manifest a parametric instability. 

Because the coordinates and velocity of a regular orbit are periodic, the frequency Q(s) 
must also be periodic, so that eq. (7) reduces to an oscillator equation with periodic 
modulation, e.g., a Hill equation. If the unperturbed geodesic is to be stable, so that it 
can exist as a regular orbit, it must be that solutions to this equation represent bounded 
oscillations, a condition that implies nontrivial restrictions on the time-dependence of Q. If 
these restrictions are not satisfied, the geodesic must instead correspond to a chaotic orbit. 

As observed by Cerruti-Sola and Pettini [7], this intuition is particularly simple to 
implement for two-dimensional systems. In this case, K a b corresponds to a symmetric 2x2 
matrix, one eigenvalue of which is necessarily zero, this corresponding to neutral stability 
with respect to infinitesimal translations along the orbit from x a (s) to x a (s + 5s). It follows 
that, in interpreting the behaviour of a small perturbation, it suffices to restrict attention 
to perturbations in the single direction orthogonal to the velocity u a , so that one is reduced 
to a single scalar equation. 

More explicitly, by exploiting the fact that, in two dimensions, the Riemann tensor has 
only one independent nonzero component, say R x yxy , it is easily seen that 

/ R x <v y ) 2 —R x v x v y \ 

\—Ry X yU U y R y X y{u ) / 

whence follows that there are two eigenvalues, namely A = and A = R x yxy [{u x ) 2 + (u y ) 2 } = 
K x x + K v y . Transforming from u a to the physical momentum p a and recalling that W is 
equal numerically to the kinetic energy, it follows that the nonzero eigenvalue 

A = Rxyxy /W 2 , (9) 

where, explicitly, 
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(10) 



The component of £ a orthogonal to u a thus satisfies 

if = -*«■ = -5*" (») 

where K = K x + K y y and R denotes the scalar curvature. Alternatively, translating back 
into physical time t one finds (cf. [7]) 

d 2 ^ 1 dWd^i 
~aW ~ W dt dt 



= -IRxyxytx = ~W 2 K^. (12) 



Cerruti-Sola and Pettini [7] have studied representative orbits in one prototypical two- 
dimensional potential, namely the Henon-Heiles potential [8], demonstrating thereby that 
one can effect a translation between various orbital properties as viewed in the ordinary 
Hamiltonian language and as viewed in this geometric language. However, it is also useful 
to study the statistical properties of ensembles of orbit segments since this facilitates a 
search for bulk regularities connecting different properties of representative orbits. Thus, in 
particular, such an investigation can provide important information about how quantities 
like short time Lyapunov exponents x(t) [9] correlate with properties of the curvature K a b 
as evaluated along some orbit. This is of interest for chaotic orbits, where one knows 
that changing values of short time Lyapunov exponents can reflect phase space transport 
through cantori [10] and other topological obstructions [11] and/or the overall complexity of 
an orbit segment, as probed by its Fourier spectrum [12]. This is also useful for regular orbits 
where, for different initial conditions, the short time exponent x(t) can converge towards 
the asymptotic Lyapunov exponent Xoo = x(t ~^ °°) at very different rates. Indeed, for 
fixed time t the value of x(t) f° r different regular orbits with the same energy E can vary 
by more than an order of magnitude. 

Equations (11) and (12) might suggest that the natural quantity upon which to focus 
is K or W 2 K. However, when realised as a function of W and its derivatives, both these 
quantities involve (cf. eq. [10]) division by positive powers of W, which, for small values 
of W, can lead to skewed statistics for a finite sampling. For this reason, it was discovered 
that cleaner results were obtained by focusing on 

K = W 3 K = WR xyxy . (13) 

Physically this combination arises if one introduces a new time coordinate r satisfying 
ds = W^ 2 dr, so that 

dr 2 2W dr dr ?1 1 ' 

The fact that the statistical properties of K = W 3 K correlate very strongly with proper- 
ties of the unperturbed orbit suggests that it is the coefficient of £i, rather than its time 
derivative, that is responsible for much of the orbit's observed behaviour. 

The work described in this paper involved examining the quantity K as evaluated along 
different orbit segments, probing in particular values of the mean (K) and the dispersion 
a^. The resulting data were used to establish trends related to these quantities, including 
how different segments - both regular and chaotic - fit into the {K)-a^ plane, and how 
the values assumed by these quantities depend on other orbital characteristics, e.g., on 
whether a regular orbit is a box or a loop or whether a chaotic segment looks nearly regular 
or particularly complex. Especially striking were correlations between (K) or and the 
values of short time Lyapunov exponents %(t) computed for the same segments. Given that 
the value of x ls strongly correlated with the overall complexity of the orbit, as probed by 
its Fourier spectrum [12], such correlations also connect statistical properties of K with the 
shape of the orbit, as viewed in configuration space. 

As noted by Pettini (private communication), the preceding justification for focusing 
K = W 3 K, rather than W 2 K, is potentially suspect since, at least in principle, any non- 
trivial time reparameterization can significantly alter the stability properties of geodesies. 
Fortunately, however, there is another interpretation which may perhaps be more easily jus- 
tified: Computing the average of W 2 K along a geodesic can be interpreted as involving a 
ratio of integrals J y/g dxdyW 2 K / J ^fg dxdy, where g denotes the determinant of the met- 
ric gij and the integration extends over the regions of the manifold along which the geodesic 



moves. However, it follows from eq. (2) that this reduces to / dxdyW 3 K / J dxdyW, and, 
to the extent that / dxdyW is approximately constant (as a consequence of virialisation) , 
one is effectively averaging the quantity W 3 K. 

The results presented below derive from an analysis of orbit segments in two different 
representative two-dimensional potentials. The first of these, 

V(x, y) = -{x 2 + y 2 ) + \{x 2 + y 2 ) 2 - \x 2 y 2 \ (15) 

corresponds to the so-called dihedral potential of Armbruster et al [13], for one particular 
set of parameter values. The second, 

V(x, y) = - (x 2 + y 2 ) + x 2 y - ^y 3 + X -x 4 + x 2 y 2 + X -y 4 

+x 4 y + \x 2 y 3 - -y 5 + \x & + x 4 y 2 + \x 2 y 4 + ^-y 6 , (16) 
3 3 5 3 45 

represents the sixth order truncation of the Toda [14] lattice potential (recall that the 

Henon-Heiles potential can be derived as the third order truncation of the Toda potential). 

These are very different qualitatively but, nevertheless, much of the observed behaviour 

is very similar for orbit segments in both potentials. In the truncated Toda potential, K 

is always non-negative. However, for certain energies in the dihedral potential K can be 

negative along parts of some orbits, although K tends to be positive most (> 75 — 80%) of 

the time. 

Ensembles of orbit segments with fixed energy E were generated by sampling the E= 
constant hypersurface and then evolving typically for a total time t = 256, this in units 
where a typical crossing time (i.e., the time required for an orbit to cross from one side 
of the potential to the other) t cr ~ 1 — 2. This is a reasonable time interval to consider 
because chaotic segments in these potentials tend to exhibit significant qualitative changes 
on a time scale ~ 100 — 200t cr [15]. However, it was verified that similar results obtain for 
somewhat longer and shorter total times. In most cases, the initial ensemble was generated 
by setting x = 0, uniformly sampling the energetically allowed portions of the y-p y plane, 
and then computing an initial p x > as a function of x, y, p y , and E. The orbits were 
integrated using a fourth order Runge-Kutta scheme with a time step St = 10 -4 . 

Strictly speaking, when extracting statistical properties of K it is most natural to anal- 
yse a time series which records relevant quantities at fixed intervals of geodesic time s, 
rather than at fixed intervals of "physical" time t, which is achieved most easily by solving 
the geodesic equation associated with W5 a ^ rather than the original Hamilton equations. 
However, it was found that the basic correlations involving quantities like (K) and a R were 
equally apparent for both sorts of time series. The discussion here focuses primarily on 
data recorded at fixed intervals 5t. This has the advantage that the conclusions derived 
here for the dihedral and truncated Toda potentials can be easily tested for other potentials, 
without the inconvenience of explicitly reformulating the evolution as a geodesic flow. 

Section II describes various trends and correlations observed in the numerical experi- 
ments, demonstrating in particular the existence of striking regularities connecting quan- 
tities like (K), and \, many of which can be interpreted in terms of other physical 
properties of the orbits. Section III summarises the principal conclusions and then shows 
that, not surprisingly, many of these can in fact be understood in terms of a simple one- 
dimensional Mathieu equation. 

II. OBSERVED CORRELATIONS AMONGST (K), a R , AND X 
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A. Correlations between (K) and 

In most, albeit not all, cases, i.e., for most energies in both potentials, the values of the 
mean (K) and the dispersion of different orbit segments are strongly correlated. Rather 
than filling a large portion of the (K)-a^ plane, the orbit segments tend to fall, at least 
approximately, along a single curve. For the dihedral potential, this curve is typically quite 
thin; for the truncated Toda potential, it can be significantly thicker. Moreover, this curve 
is typically characterised overall by a negative slope, so that orbits with larger (K) have 
smaller a^. 

Figures 1 and 2 exhibit plots of (7^ as a function of (K) for several different energies 
in, respectively, the dihedral and truncated Toda potentials. Figures Id - h and Figs. 2c 
and d, each characterised by a single curve, perhaps with one or more intermittent gaps 
(as in Figs, le and f) or a slight "wiggle" (as in Fig. lh) represent typical behaviour. 
The fact, manifest visually, that the values of (K) and <r^> are strongly correlated can be 
quantified by computing the rank correlation 1Z((K), — a^). For these typical cases in the 
dihedral potential, the correlation between (K) and — is usually very strong, 1Z > 0.98 
or more. For the truncated Toda potential, the correlation is often somewhat weaker, 
K((K), -a k ) ~ 0.9 - 0.95, but still significant. 

This sort of correlation between (K) and a p., seemingly suggestive of ordered behaviour, 
is perhaps not surprising for regular orbits, where the motion is multiply periodic. Thus, 
e.g., it is easy to envision a sequence of regular box or loop orbits characterised by smoothly 
varying values of (K) or o^. However, such correlations might seem less expected for chaotic 
orbits where the motion is aperiodic. It is therefore significant that, when viewed in such a 
(K) -a^ plot, regular orbits need not stand out in any obvious way. Consider, e.g., the orbit 
ensemble with E = 4.0 used to generate Fig. If. Here one finds that there are three distinct 
types of regular orbits, a large number of orbits with 59 < (K) < 63 and 9 < < 26, 
a large number of orbits with 54 < (K) < 57 and 33 < < 36, and a small number of 
orbits with 43 < (K) < 45 and 49 < < 50. It is clear that all the points below the gap 
in Fig. If are associated with regular orbits, and a careful examination of the data points 
allows one to distinguish minute differences between the locations of the regular and chaotic 
orbits in the intermediate regime, 54 < (K) < 57 and 33 < < 36. However, it is evident 
that, overall, the regular and chaotic orbits coexist along a relatively narrow curve [16]. 

This can be associated tentatively with the fact that, even though chaotic orbits are 
intrinsically aperiodic, finite chaotic segments often manifest a fair amount of regularity. 
Thus, e.g., as discussed more carefully below, one observes oftentimes that the power spec- 
tra, \x(oj)\ and for a chaotic segment typically appear visually to be constructed 
from "pieces" appropriate for a small number of regular orbits [12,17]. This interpretation 
is especially natural given the recognition that the most striking exception to the simple 
pattern described hitherto is associated with very low energies in the dihedral potential, 
where few, if any, regular orbits exist. 

Figure la, a seemingly structureless set of points completely different from the lines 
observed in the remaining panels of Figs. 1 and 2, derives from an ensemble of segments 
with E = —0.05, an energy where a sampling of nearly 1000 different initial conditions 
yielded only chaotic orbits. At slightly higher energies, regular orbits begin to appear and, 
unlike the chaotic orbits, they seem concentrated largely along lines in the (K)-a^ plane. 
Thus, e.g., as illustrated in Fig. lb, for E = +0.05 there are four different types of regular 
orbits, three of these associated with the three conspicuous lines and the fourth associated 
with the largest values of above the central chaotic region. As the energy is raised to yet 
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higher values, the regular families associated with the two upper lines eventually disappear, 
so that, e.g., as illustrated in Fig. Id for E = 1.0, all the orbits fit approximately onto a 
single curve, with the chaotic orbits at high and the regular orbits at low a^. When 
the energy is raised to a higher value E ~ 3, other regular orbit families appear, leading 
eventually to the aforementioned behaviour at E = 4.0. 

The observed behaviour at low energies in the truncated Toda potential is quite different, 
this presumably reflecting the fact that, in this case, at very low energies there is no global 
stochasticity. For energies E < 5.0 or so, all the orbits present, both regular and chaotic, fit 
into two distinct lines in the (K) -a^- plane. Consider, e.g., the energy E = 0.5 exhibited 
in Fig. la. Here most of the upper line is occupied by loop orbits which manifest a discrete 
27r/3 rotation symmetry, whereas the lower line is occupied completely by box and banana 
orbits that break this symmetry. There are only a very few chaotic orbits at this energy, 
and all of them fit at the high end of the upper line. As the energy increases, the lower 
line eventually shrinks and merges into the upper line, so that, ultimately, different families 
of regular orbits coexist with chaotic orbits along a single thickened curve. 



Ordinary Lyapunov exponents x> which [18] probe the average instability of some trajec- 
tory in an asymptotic t — > oo limit, manifest a fundamental distinction between regular and 
chaotic orbits. For regular orbits all the Lyapunov exponents vanish identically, whereas 
chaotic orbits have at least one x that is positive. Orbits in a D-dimensional system have 
2D distinct Lyapunov exponents, these corresponding to perturbations in 2D independent 
phase space directions. For a time-independent Hamiltonian system, two of these exponents 
must vanish (this reflecting neutral stability with respect to perturbations that translate an 
orbit from x a (t) to x a (t+5t) and to perturbations orthogonal to the constant energy surface) 
and the remaining exponents must come in pairs, ±%. It follows that, for two-dimensional 
Hamiltonian systems, the only fundamental distinction is between regular orbits, for which 
all the x's vanish, and chaotic orbits, which have one exponent x > 0. 

However, one can also introduce short time Lyapunov exponents x(t) which provide 
information about the average instability of orbit segments over a finite interval. Thus, in 
particular, for any phase space perturbation Sz, one can define [9] 



where |.| represents a suitable norm. For a generic initial perturbation, this x{t) will 
converge towards the largest Lyapunov exponent in the limit t — > oo, independent of the 
detailed choice of norm. By contrast, at finite times the computed will depend on 
both the initial perturbation and the choice of norm. Suppose, however, that | . | is taken 
as the natural L 2 phase space norm, i.e., ||<5Z|| 2 = (Sx) 2 + (Sy) 2 + (5p x ) 2 + (5p y ) 2 . In this 
case, one knows that the computed x(i) will be insensitive to the detailed choice of initial 
perturbation for times t 3> l/x(i). 

It follows that, for chaotic segments integrated in the dihedral and truncated Toda 
potentials for times as long as t = 256, the computed x(t) is nearly independent of the 
initial 5Z. However, as described below the values of x(t) computed for regular orbits 
can, and in certain cases do, exhibit a significant dependence on bZ . For the experiments 
described in this paper, %(t) was computed [18] by introducing a small perturbation of 
magnitude \ \6Z\ \ = 10 -10 , evolving both unperturbed and perturbed initial conditions, and 



B. Correlations between curvature and %(£) 




\\8Z(t)\\ 
\\SZ(0)\\ 



(17) 



7 



periodically renormalising the perturbation to an amplitude \\SZ\\ = 10~ 10 at intervals 
At = 10. Unless stated otherwise, the initial perturbation was taken as SZ = 5x = 10~ 10 . 

The objective here is to show that, both for regular and chaotic orbit segments, strong 
correlations exist between the value of the short time and such properties of the curva- 
ture as (K) and a^. Because the quantities (K) and are themselves correlated, it would 
seem equally reasonable to look for correlations between and either (K) or (or any 
combination of these two quantities). For specificity, most of the discussion will focus on 
correlations between and x(t), although several Figures exhibit examples of correlations 
between (K) and x(i). Figures 3 and 4 exhibit plots of % as a function of cr^-, generated 
respectively for the dihedral and truncated Toda potentials for the same ensembles used to 
generate Figs. 1 and 2. 

Consider first the case of orbit ensembles evolved in the dihedral potential. For E = 
—0.05, where there are few if any regular orbits, a plot of x as a function of shows 
little obvious structure: all that one sees is a seemingly random scattering of points at 
values of x well separated from x = 0. However, as E increases, one begins to observe the 
existence of regular regions, these corresponding to ranges of where, for some segments, 
X assumes values much smaller than the values associated with chaotic segments. Thus, 
e.g., as illustrated in Fig. 3b, for E = 0.05 one observes an extended low band that 
corresponds to one regular line in Fig. lb, a collection of points near = 6.0 corresponding 
to the second line, another collection near = 7.4 corresponding to the third line, and 
a small number of small x points near = 7.8 corresponding to the high points in 
Fig. lb. 

For E > 0.2 or so, a simpler pattern emerges which includes only two types of regular 
orbits, these concentrated at especially low and high values of (or, equivalently, high 
and low values of (K)). At somewhat higher energies, the high family disappears, only 
to be replaced by a new regular family concentrated at intermediate values of a^. At yet 
higher energies, one sees two large regular regions, concentrated at the largest and smallest 
values of a^, along with some intermediate regions with small x's, these associated with 
nearly vertical lines in the cr^-X plane. Viewed in a surface of section, the low regular 
region corresponds to a large island of loop orbits; the high region corresponds to a 
large island of box orbits. The regular orbits associated with the vertical lines correspond 
to smaller islands embedded in the stochastic sea. This structure appears to persist to very 
high energies. 

For these relatively high energies, E > 6.0 or so, a plot of x as a function of or (K) 
exhibits several striking regularities. One significant point, well illustrated for E = 10.0 
in Fig. 5, is that the transition from regularity to chaos observed at low and high values 
of Ofr (or high and low values of (K)) is relatively abrupt. Thus, e.g., the small a^, 
large (K) loop orbits can be viewed as a sequence beginning at « 22 that terminates 
at a value ~ 112, whereas the large a^, small {K) boxes can be viewed as a sequence 
extending upwards from (K) ~ 124 and terminating at a value (K) ~ 135. The chaotic 
segments situated near the boundary with the outer regular regions tend typically to be 
"confined" or "sticky" chaotic orbits [19] trapped near the regular regions by cantori [10] 
which, oftentimes, only escape to travel throughout the remainder of the stochastic sea on 
a time scale t > 256. It is also apparent from Fig. 5 that (even away from the boundaries) 
chaotic segments with values of furthest from the high and low a ^ regular regions tend 
systematically to have larger %'s than do chaotic segments with values of closer to the 
outer regular regions. The net result is a curve which, viewed broadly, resembles a Greek 
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A. 

The overall shape of the ap, — % curve in the regular region depends on the choice of 
the initial seed SZ used in computing the short time x(*)- If> as i n Figs- 5a and b, one 
computes % from a seed SZ = Sx, it is apparent that the segments with values of (or 
(K)) closest to the central chaotic region, 112 < < 142 or so, tend overall to have larger 
values of x than do orbits with values of further from this chaotic region. However, if 
X is generated instead from a seed with a significant component in one of the other three 
phase space directions, this trend is significantly diminished. This is illustrated in Figs. 5c 
and d, which exhibit generated for the same initial conditions from a seed SZ = Sy. 
This difference presumably reflects the fact that, since the initial conditions were sampled 
from an x = surface of section, a perturbation SZ = Sx tends to be more nearly aligned 
along a direction of neutral stability, namely translation from x(t) to x(t + St), than do 
perturbations with a nonzero Sy. 

At lower energies, where the high region is absent or not well developed, the right 
side of the A is missing. However, one still observes that chaotic segments with closer 
to the regular region tend to have smaller values of x than do chaotic segments far from the 
regular regions; and, at least for a seed SZ « Sx, that regular orbits closer to the chaotic 
region tend to have larger values of x- 

One obvious complication associated with the common A pattern, well illustrated in Fig. 
5, is the presence of one or more nearly vertical lines in the chaotic region, extending from 
very low to relatively high values of x- The lowest values of x seem too small to be associated 
with chaotic segments, but the upper values seem too large to be associated with regular 
orbits. In fact, these lines contain two different classes of orbits, namely regular orbits, 
for which %(t) eventually decays to zero, and confined chaotic orbits, trapped temporarily 
near a small regular island, which eventually escape through one or more cantori to travel 
unimpeded through the stochastic sea. Viewed in configuration space, the regular orbits 
corresponding to these lines correspond to periodic orbits confined to an annulus or a "figure- 
eight-shaped" region. The confined chaotic orbits correspond to aperiodic orbits which, for 
a long time, are trapped in almost the same region (whence follows the fact that they have 
nearly the same (K) and as do the regular orbits), but eventually escape to probe the 
remaining chaotic phase space regions. 

Evidence for these assertions is provided in Fig. 6, which summarises an investigation 
of the longer time evolution of the initial conditions which led to the near vertical line in 
Fig. 5b at ~ 139. A random sampling of 1261 initial conditions with E = 10.0 evolved 
for a total time At = 256 led to 81 segments with 138 < a^- < 140 and x(At) < 0.1. These 
81 initial conditions were subsequently integrated for a significantly longer time, t = 8 x 256, 
with x(t) recorded at regular intervals, and the resulting orbits partitioned into 8 segments 
of length At = 256, for which values and (K) were computed. The recorded values of 
x(t) were then analysed to extract short time Lyapunov exponents for successive intervals 
of length At = 256, the exponent for the segment extending from t k = kAt to t k + At being 
defined via the obvious relation (cf. [12]) 

(A . x (tk + At) x (t k + At) - t kX {t k ) ns , 
X\At k ) = — . (18) 

The tiny dots in all four panels of Fig. 6 exhibit the values of x an d cr^ generated for 
the original ensemble of 1261 orbits evolved for a time t = 256. The diamonds in Fig. 6a 
highlight the 81 segments concentrated initially along the line near = 139. The diamonds 
in Fig. 6b show the values of x an d cr^ derived from the same 81 initial conditions for the 
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period 256 < t < 512. Figures 6c and d extend the results to later intervals 512 < t < 768 
and 768 < t < 1024. It is clear that, as time elapses, the larger x diamonds escape from 
the line and move to other portions of the chaotic regions, whereas the smaller \ diamonds 
evolve closer to % = 0. 

The other complication common to the A pattern is the existence of smaller scale struc- 
tures in the regular regions. As noted already, for short time exponents generated from a 
seed 5Z directed nearly in the x-direction, many/most regular segments fit along a curve 
with x decreasing as one moves away from the central chaotic region to much larger or 
smaller values of a p.. Alternatively, for more generic choices of SZ, a plot of % as a function 
of Gj^ exhibits an upper envelop that is is more nearly independent of <7j>. However, in each 
case this upper curve is not the whole story. 

As illustrated in Figs. 5c and d, for generic values of 5z the low loop orbits exhibit 
two additional features, namely (1) small scale oscillations in the cr^-x or (K)-X plane, 
and (2) an excess probability of finding segments with especially low values of x closer to, 
rather than further from, the central chaotic region. Alternatively, as illustrated in Figs. 5a 
and b, for Sz oriented more nearly in the x-direction one can identify various "sub-families" 
of regular orbits which, for a fixed value of cr^-, tend to have much smaller values of % than 
other regular orbits with (nearly) the same value of a^. Thus, e.g., if one focuses on a fixed 
interval of values of and compares regular orbits with larger values of x that fall along 
the upper line with regular orbits with smaller values of x, ne or she finds typically that 
there is a relatively clear cut difference between the two sets of orbits, as probed by the 
detailed shapes of their power spectra or, in some cases, by the overall shape of the orbit. In 
many cases, the regular orbits that have smaller values of x tend to be "simpler" or "more 
regular" in appearance. This is illustrated in Fig. 7, which exhibits two representative 
regular orbits extracted from the interval 51.0 < < 58.0. Figure 7a shows a typical 
low-x orbit, whereas Fig. 7b shows a typical orbit lying along the upper curve. It is obvious 
that the higher- x orbit is (nearly) space filling within an annulus, whereas the lower orbit 
exhibits more structure. 

The detailed patterns relating x to ff^ are different for initial conditions evolved in 
the truncated Toda potential, but the data still admit to a similar interpretation. At 
relative low energies, E ~ 2 — 8, one sees patterns qualitatively similar to those observed 
for E > 6 or so in the dihedral potential. Alternatively, for somewhat higher energies, 
E > 25 or so, one sees instead a pattern reminiscent of those observed at relatively low 
energies in the dihedral potential where, in the absence of a significant population of large 
o~x box orbits, one observes simply that, for chaotic segments, larger correlates with 
larger x- Superficially, the behaviour for E = 20.0, illustrated in Fig. 4c, might seem quite 
different from what is observed in the dihedral potential since here the chaotic segments 
with the smallest values of x are those furthest from the regular regions. However, when 
placed in an appropriate context this behaviour is not surprising. At energies slightly above 
E = 20.0, this low x chaotic region merges into a relatively large regular island, which is 
well established by E = 25.0. For E = 30.0, one sees three large islands at low a^, the 
two corresponding to the lowest values of very similar qualitatively to the islands for 
E = 10.0 in the dihedral potential associated with the two vertical lines. 

C. Curvature and Complexity 

The preceeding shows that there is a strong correlation between the overall stability or 
instability of an orbit segment, as probed by the value of a short time Lyapunov exponent 
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x(t), and the values of (K) and associated with that segment. However, earlier work [12] 
has shown that, for chaotic segments, there also exists a strong, nearly linear, correlation 
between %(£) and the complexity of the segment, as probed by the form of its Fourier 
spectrum. One thus anticipates that there should exist correlations between (K) or and 
the complexity of the segment. 

As discussed more carefully in Ref. [12], n(k), the complexity of an orbit segment at 
threshold k, can be defined in the following way: Given the values of the configuration 
space coordinates, say x and y, at fixed time intervals At, sampled, e.g., ~ 5 — 10 times per 
crossing time, compute the power spectra, \x(u)\ and \y(uj)\ in the usual way [20]. Next 
determine the minimum number of frequencies, n x (k) and n y (k), required respectively to 
capture a fixed fraction k of the x- and y-powers. The total complexity is then defined as 

n(k) =n x (k) +n y (k). (19) 

This notion of complexity is motivated by the idea that "comparatively regular" seg- 
ments have most of their power concentrated at or near a few special frequencies, whereas 
"wildly chaotic" segments have spectra with significantly broader band power. Experience 
with ensembles of orbits evolved in several different potentials has indicated that values 
k ~ 0.9 — 0.95 yield relative complexities in good agreement with subjective impressions 
based on visual inspection of segments in the x — y plane. The results described in this 
paper were derived from integrations where coordinates were recorded at fixed intervals 
At = 0.125 for a total time t = 256, this leading to 2048 points and hence a Fourier series 
with 2048 frequencies. 

As stated already, for orbit segments in both the dihedral and truncated Toda potentials, 
as well as for other systems [12], there is a strong, nearly linear correlation between %(t) 
and quantities like n(0.9) or n(0.95). For ensembles where x an d n(k) both assume a broad 
range of values, the rank correlation typically assumes a value lZ(x,n(k)) > 0.85 — 0.9. If, 
alternatively, most of the segments are concentrated in the same part of the x ~ n (k) plane, 
1Z can be smaller but still remains appreciable. 

This correlation can be understood intuitively in terms of the fact that chaotic orbit 
segments often appear visually to be comprised of "pieces" of various regular orbits that 
exist at the same energy [16]. Consider, e.g., orbits with E = 10.0 in the dihedral potential. 
Here one finds that chaotic segments with values of near those appropriate for the two 
large regular regions tend overall to look more regular, and to have smaller complexities, 
than segments located at values of further from the regular regions. Moreover, chaotic 
segments near the low regular region tend to look "loopy" and to have power spectra 
dominated by features appropriate for the loop orbits that exist at small <jj>; and similarly, 
chaotic segments near the high regular region tend to look "boxy" and to have spectra 
similar to the high ap, box orbits. 

The degree to which a quantity like n(0.9) correlates with x, both in itself and in relation 
to other quantities like (K) or a^, can be gauged from Fig. 8, which exhibits data generated 
in the dihedral potential for E = 10.0 and E = 1.0. The two left panels exhibit plots of 
X as functions of (K), whereas the two right panels plot n(0.9) as a function (K). It is 
clear that, at least in the chaotic regions, the plots of x an d n(0.9) exhibit the same basic 
features, although the plots of n(0.9) seem more "blurred." 

There is, however, a significant difference for the regular orbits. Plots of n(0.9) (or 
n(0.95)) as a function of or (K) exhibit large systematic oscillations not manifested in 
the plots of x, whereas the small scale structures present in the plots of x (cf. Fig. 5) arc 
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absent in plots of n(0.9). Consider, e.g., E = 10.0. Here all the segments at large (K) are 
regular loops, but the loops near (K) = 175 are clearly special in that n(0.9) assumes values 
that are particularly small. This reflects the fact that, in this region, peaks in the spectra 
\x(uo)\ and \y(u>)\ are much sharper than at somewhat larger and smaller values of (K). 

This is a finite sampling effect, reflecting the fact that the orbits segments were inte- 
grated for a relatively short time t = 256. If the total time t is changed, but data are still 
recorded and analysed at 2048 equally spaced intervals, one finds that, as would be expected, 
the values of (K) and computed for regular segments remain essentially unchanged. By 
contrast, however, the details of the oscillations can be alterred significantly. Thus, e.g., for 
E = 1.0 the location of the dip near (K) = 17.5 will move and other minima can appear. 
Moreover, one observes that, as t increases, the computed complexities decrease and the 
amplitudes of the oscillations damp. 

One final correlation remains to be stated. To say that a regular orbit is close to the 
chaotic region is a statement both about its location (say) in the cr^-x plane and its actual 
location in configuration space, as probed, e.g., by a surface of section. For example, orbits 
in the large regular regions with values of (K) and a p. further from the central chaotic 
region tend systematically to lie closer to the center of a regular island than do regular 
segments with values of (K) and closer to the central chaotic region. Similarly, the 
vertical lines associated with regular and confined chaotic orbits are associated with smaller 
regular islands embedded in the stochastic sea. 

These facts are illustrated in Fig. 9. Here the small dots represent successive sections 
generated from a single "unconfined" chaotic initial condition that systematically avoids the 
phase space regions near the two large regular islands. The remaining symbols represent 
different initial conditions corresponding to regular and "confined" chaotic orbits generated 
from a uniform sampling of the y-p y plane. The smaller of the two large islands, located 
near y = p y = and populated with stars, corresponds to box orbits. The larger, displaced 
from y = 0, corresponds instead to loop orbits. Here the initial conditions were binned 
into three classes. The diamonds represent segments with a p. assuming particularly small 
values, whereas the squares represent values of especially close to the chaotic region. 
The triangles correspond to intermediate values. The remaining pluses and crosses represent 
initial conditions which led to the two nearly vertical lines exhibited in Fig. 5. 

III. SUMMARY AND INTERPRETATION 

Several general trends emerge from the experiments described in the preceeding Section. 
Most striking and fundamental, perhaps, is the fact that, for an ensemble of orbit segments 
of fixed energy E, a plot of the dispersion function of the mean (K) usually 

assumes a very simple form. In general, regular and chaotic segments appear to coexist 
along a single curve in the (K)-ap- plane, although for both potentials one sees more 
complicated structures at very low energies. The fact that chaotic segments fit along a 
single curve indicates that, even though they are aperiodic, they still exhibit significant 
statistical regularities. 

Plots of the short time Lyapunov exponent x(t) as a function of (K) or up, also exhibit 
simple, distinctive patterns. For example, one finds that, largely independent of the initial 
perturbation 5Z used in calculating x> chaotic segments with values of (K) (or a p.) far 
from the regions associated with regular orbits tend to have larger values of x(t) than 
do chaotic segments with values of (K) (or a^) nearer the regular regions. For regular 
orbits, the computed values of x(*)> an d hence any correlations with (K) or a p., depend 
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more sensitively on the initial perturbation. For most initial perturbations, it seems that, 
overall, the typical x(t) is nearly independent of (K) and a p.. However, for some classes of 
perturbations one finds instead that segments with values of (K) (or a^) further from the 
chaotic regions tend to have smaller values of x than do regular segments with values closer 
to the chaotic region. 

The correlations between x(t) and {K) (or a^) observed for chaotic orbits are very 
similar to the observed correlations between the complexity n(k) and (K) (or o^)- This 
reflects the fact that chaotic segments with values of (K) and closer to the regular 
regions tend to look more regular and, consequently, to have less complex Fourier spectra 
than do orbits with values further from the regular region. 

One other common feature, observed for both complexities and short time exponents, is 
the presence of one or more nearly vertical lines in (say) the - x or - n(k) plane. These 
are comprised of two distinct types of segments, namely segments of regular orbits trapped 
forever in an island of stability by invariant K AM tori and "confined" chaotic orbits, trapped 
temporarily by cantori in a similar configuration space region which, however, eventually 
diffuse through the cantori to probe the remaining portions of the chaotic sea. 

Many of these qualitative features can be reproduced in a much simpler context, where 
the linearised scalar Jacobi equation (11) is replaced by a linear Matthieu equation of the 
form [21] 

^ = -(u, 2 + 7 sini)6, (20) 

where u> 2 and 7 play the roles, respectively, of (K) and a^. 

If one studies solutions to eq. (20) as a function of 0J 2 and 7 (or, by analogy, (K) and 
afr) he or she finds that the u - 7 plane is partitioned into a large number of disjoint regions, 
corresponding respectively to "regular" and "chaotic" orbits. In the regular regions, it is 
apparent that, in the t — > 00 limit, the analogue of the ordinary short time Lyapunov 
exponent, 



X(t) = \ In 



\\SZ(0)\\ 



(21) 



with ||5^|| 2 = |£i| 2 + |£i| 2 , eventually converges towards zero. Alternatively, in the chaotic 
regions x(t) converges towards a time-independent positive value. It follows that a generic 
curve passing through the cj-7 (or (K)-a^) plane, analogous to the curves exhibited in 
Figs. 1 and 2, will include both regular and chaotic regions with sharp transitions between 
the two. 

For values of oj and 7 far from the chaotic portions of the curve, one finds typically 
that, for reasonably large t, the computed value of x is relatively insensitive to the detailed 
perturbation and to the exact values of to and 7. However, for values of uj and 7 closer to 
the chaotic regions, more complicated patterns can arise. Viewed in an asymptotic t — ► 00 
limit, the transition from regular to chaotic is abrupt. However, the transition is smooth 
in the sense that, when evaluated for some finite time t, the short time x(t) computed for 
"chaotic" values of u and 7 close to the regular regions is typically smaller than the %(t) 
computed for values that are further from the regular region. This is consistent with the 
fact, manifested in Figs. 3-5, that orbits with values of (K) and further from the 
regular regions tend to have larger short time Lyapunov exponents. Given this observation, 
it is also easy to understand (cf. Fig. 6) why, for "sticky" chaotic segments confined near 
regular regions by cantori, the computed %(t) tends to be smaller than for other chaotic 
segments that travel further from the regular regions. 
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A detailed discussion of the sense in which the onset of chaos can be understood in 
terms of a more general Hill equation has been presented by Cerruti-Sola and Pettini [7] in 
the context of their analysis of individual orbits. 

The work described in this paper has established the existence of a strong correlation 
between such geometric properties of chaotic orbit segments as (K) and ap, and the short 
time Lyapunov exponent \. However, one might expect that this local, short time analyses 
could also be extended to establish a global connection between moments of K, as defined 
along an infinite chaotic geodesic, and the ordinary Lyapunov exponent Xoo, as defined in a 
t — > oo limit. Casetti et al [22] have shown that, at least for systems with a large number N 
of degrees of freedom, such a connection does indeed exist. Specifically, by assuming that 
the curvature experienced by a chaotic orbit can be approximated as a random process, 
characterised by a mean (K) and a dispersion ax, they obtained an analytic approximation 
to the largest Lyapunov exponent which, in at least some cases, agrees extremely well with 
numerical computations. 

It would seem unlikely that such a simple Gaussian approximation, which these authors 
motivate in the large N limit, will work well for the very special case N = 2. However, one 
might nevertheless hope that, in some fashion, the value of the positive Lyapunov exponent 
for a two-dimensional system can again be related to the moments of K or K evaluated 
along the chaotic orbit. This possibility is currently under investigation. 
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FIGURE CAPTIONS 



FIG. 1 - The mean (K) and dispersion for orbit ensembles with eight different energies 
E evolved in the dihedral potential for a total time t = 256. (a) E = —0.05. (b) E = 0.05. 
(c) E = 0.25. (d) E = 1.0. (e) E = 2.0. (f) E = 4.0. (g) E = 8.0. (h) E = 20.0. 

FIG. 2 - The mean (K) and dispersion a p. for orbit ensembles with four different energies E 
evolved in the truncated Toda potential for a total time t = 256. (a) E = 0.5. (b) E = 3.0. 

(c) E = 20.0. (d) E = 30.0 

FIG. 3 - The short time Lyapunov exponent \ plotted as a function of ap~ for the same 
orbit ensembles as in Fig. 1. 

FIG. 4 - The short time Lyapunov exponent x plotted as a function of up- for the same 
orbit ensembles as in Fig. 2. 

FIG 5 - (a) Short time Lyapunov exponents x(i), generated with an initial perturbation 
Sx = 10~ 10 and plotted as functions of op- for an ensemble of orbit segments with E = 10.0 
evolved in the dihedral potential for a total time t = 256. Only values of \ < 0.16 are 
shown, (b) The same xW s plotted as a function of {K). (c) The analogue of FIG. 5a, 
generated for the same ensemble, but now allowing for an initial perturbation Sy = 10~ 10 . 

(d) The analogue of FIG. 5b, generated for an initial perturbation Sy = 10~ 10 . 

FIG. 6 - The evolution of regular and confined chaotic orbits for E = 10.0 in the dihedral 
potential, with \ plotted as a function of a^. The dots represent analogues of FIG. 3, with a 
total integration time t = 256. The diamonds represent the same initial conditions analysed 
over different time intervals, (a) < t < 256. (b) 256 < t < 512. (c) 512 < t < 768. (d) 
768 < t < 1024. 

FIG. 7 - Two orbits with E = 10.0 evolved in the dihedral potential for a time t = 256. 
(a) An orbit with (K) = 183.4, a K = 56.4, Q = 0.308, and X = 0.0051. (b) An orbit with 
(K) = 181.6, g k = 62.7, Q = 0.345, and \ = 0.0116. 

FIG. 8 - (a) The short time exponent \ plotted as a function of (K) for a collection of orbits 
with E = 10.0 evolved in the dihedral potential for a total time t = 256. (b) The complexity 
n(0.9) plotted as a function of (K) for the same set of segments, (c) The analogue of FIG. 
8 (a) for a collection of segments with E = 1.0. (d) The analogue of FIG. 8 (b) for the 
segments of FIG. 8 (c). 

FIG. 9 - An x-p x surface of section for E = 10.0 in the dihedral potential. The dots derive 
from a single unconfined chaotic orbit evolved for 2000 intersections. The remaining points 
represent some of the initial conditions used to generate FIGS. 5 and 6. Diamonds represent 
loop orbits with < 72.0. Triangles represent loop orbits with 72.0 < < 90.0. Squares 
represent loop orbits with 90.0 < < 110.0. Stars represent box orbits with > 141.0. 
Pluses represent points along the near- vertical line with 120.0 < ap, < 122 and x < 0.08. 
Crosses represent points along the line with 138.0 < cr^> < 140.0 and x < 0.1. 
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